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Abstract 

Motivated by constraint-based CAD software, we develop the foun- 
dation for the rigidity theory of a very general model: the body-and- 
cad structure, composed of rigid bodies in 3D constrained by pairwise 
coincidence, angular and distance constraints. We identify 21 relevant 
geometric constraints and develop the corresponding infinitesimal rigidity 
theory for these structures. The classical body-and-bar rigidity model can 
be viewed as a body-and-cad structure that uses only one constraint from 
this new class. 

As a consequence, we identify a new, necessary, but not sufficient, 
counting condition for minimal rigidity of body-and-cad structures: nested 
sparsity. This is a slight generalization of the well-known sparsity condi- 
tion of Maxwell. 
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1 Introduction 



This paper initiates the study of and sets up the foundation for the rigidity 
theory of a large class of 3D geometric constraint systems. These systems are 
composed of rigid bodies with specific coincidence, angular and distance con- 
straints and are called body-and-cad structures. To the best of our knowledge, 
these constraints have not been systematically studied before from this perspec- 
tive. 

Motivation. Popular computer aided design (CAD) software applications 
based on geometric constraint solvers allow users to design complex 3D sys- 
tems by placing geometric constraints among sets of rigid body building blocks. 
The constraints are specified by identifying geometric elements (points, lines, 
planes, or splines) on participating rigid bodies. Detecting when a user has 
created a fully-defined sub-system or has added a redundant (or inconsistent) 
constraint are important problems for providing informative feedback. How- 
ever, analyzing all constraints simultaneously is a very difficult problem. In 
this paper, we focus on a subset of these constraints that are amenable to a 
rigidity-theoretical investigation. 

Underlying classical rigidity theory results is a general proof pattern, span- 
ning algebraic geometry (for rigidity), linear algebra (for infinitesimal rigidity) 
and graph theory (for combinatorial rigidity). The ultimate goal is a full com- 
binatorial characterization of generically minimally rigid structures, but such 
results are extremely rare: 3D bar-and-joint rigidity remains a conspicuously 
open problem [7], while the 2D version is fully understood [12] . An important 
step along the way is identifying a pattern in the rigidity matrix developed as 
part of the infinitesimal rigidity theory for the structures. While this is straight- 
forward for the well-known bar-and-joint model, it is more complicated in the 
body-and-bar model. In this paper, we formulate the even more involved rigidity 
matrix for the body- and- cad model. 

Results. We define a body- and- cad structure to be composed of rigid bodies 
connected by pairwise coincidence, angular (parallel, perpendicular, or arbitrary 
fixed angular) and distance constraints. The constraints occur between specified 
points, lines or planes (called geometric elements). Besides the well-studied 
distance constraint between points (as in body-and-bar structures), we identify 
20 new pairwise constraints. We label constraints by the geometric elements 
involved, e.g., a line-plane perpendicular constraint between bodies A and B 
indicates that a line on A is perpendicular to a plane on B. The complete set of 
body-and-cad constraints that we study is further subdivided into six categories: 

• Point-point constraints: coincidence, distance. 

• Point-line constraints: coincidence, distance. 

• Point-plane constraints: coincidence, distance. 

• Line-line constraints: parallel, perpendicular, fixed angular, coinci- 
dence, distance. 
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• Line-plane constraints: parallel, perpendicular, fixed angular, coinci- 
dence, distance. 

• Plane-plane constraints: parallel, perpendicular, fixed angular, coinci- 
dence, distance. 

We develop the pattern of the rigidity matrix and identify a necessary 
combinatorial counting property called nested sparsity, which is the counterpart 
of the well-known Maxwell condition [18] for fixed length rigidity. We also show 
that this condition is not sufficient. However, it can be used as a filter for 
finding candidate rigid components. Finally, we present an efficient algorithm 
for nested sparsity, based on pebble game algorithms previously developed for 
sparse graphs. 

Related work. Classical rigidity theory [7] focuses on distance constraints 
between points [12] or rigid bodies [30] [34]. Direction constraints (where 2 
points are required to define a fixed direction, with respect to a global coordinate 
system) are well-understood and arise from parallel redrawing applications [35] . 
Motivated by CAD systems, Servatius and Whiteley present a characterization, 
which can be viewed as a generalized Laman counting property, for 2D systems 
with both length and direction constraints [24 . 

Work on angular constraints has also focused on combinatorial characteri- 
zation results. Zhou and Sitharam [37 characterize a large class of 2D angular 
constraint systems along with a set of combinatorial construction rules that 
maintain generic independence. Saliola and Whiteley [21] prove that, even in 
the plane, the complexity of determining the independence of a set of circle 
intersection angles is the same as that of generic bar-and-joint rigidity in 3D. 
A full characterization for angular constraints of the nature that appear in this 
paper is further described in [16] [13] . 

Combinatorial sparsity conditions are intimately tied with rigidity theory, 
appearing often as necessary conditions (as for bar-and-joint rigidity) and some- 
times even as complete characterizations (as for 2D bar-and-joint and body-and- 
bar frameworks in arbitrary dimension) [36j [12j [30] . Pebble game algorithms 
have been developed for solving sparsity problems [14, 29, 15] . These algorithms 
do not apply, however, to the so-called (3, 6)-counting conditions known to be 
a necessary, but not sufficient, condition for 3D bar-and-joint rigidity. In fact, 
no efficient algorithm is known for these counts. 

Related work on the constraints studied in this paper has appeared in the 
CAD research community, usually within the context of decomposition ap- 
proaches; a survey may be found in [25] [11]. In this setting, a geometric con- 
straint system (GCS) is formulated as an algebraic system of equations. Due 
to the complexity of solving such a system, it is traditionally decomposed into 
structured sub-systems that can be solved and later recombined to obtain a solu- 
tion to the original GCS. In the process of decomposition, approximate notions 
of combinatorial rigidity have been used [TO] [9] . 

Results in the CAD literature have observed that angular constraints exhibit 
special behavior. For the so-called generalized Stewart platform, [5] gives ex- 
plicit equations that highlight this distinction. Gao et al. [6 present a method 
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for analyzing 2D and 3D systems with a restricted set of coincidence, angular 
and distance constraints. Both [5 and [6] treat angular constraints separately, 
implicitly using natural necessary counting conditions to do so. We consider 
analogous systems from the rigidity theory perspective, expressing them in- 
finitesimally using Grassmann-Cayley algebra; the shape of the rigidity matrix 
described in Section [3] explicitly reveals the distinct treatment of angular con- 
straints. Grassmann-Cayley, Clifford algebras and geometric algebras often ap- 
pears in the context of CAD or geometric theorem proving; see, e.g., [T7] 123]. 
Recent work of [20] expresses constraints from a similar perspective when pro- 
viding a foundation for software to build a GCS. 

Incidence constraints have been studied previously in connection with Ge- 
ometric Theorem Proving [17] [19] for projective incidence theorems. Sitharam 
et al. [26] [27] [28] formalize the question of obtaining a well- formed and optimal 
system of algebraic equations to resolve a collection of incident rigid bodies. [26] 
studies "well-formedness," a condition necessary to avoid dependent equations, 
and a new, underlying matroid whose independent sets capture this. A combi- 
natorial measure of algebraic complexity of the system of equations is described 
in [27] , and another underlying matroid is used to optimize this measure. In [28] , 
it is shown how to reconcile the independent sets of the prvious two matroids 
to obtain an optimal, well-formed system. 

Structure. Section [2] gives a brief overview of the required mathematical back- 
ground. Section [3] develops the foundations for the infinitesimal rigidity theory, 
providing the basic building blocks used for each new constraint. Each of the 
full set of constraints is then expressed using these building blocks in Section [4] 
resulting in the complete derivation of the rigidity matrix. Section [5] identifies a 
new combinatorial property resulting from the structure of the rigidity matrix; 
this nested sparsity condition, while necessary, is shown not to be sufficient with 
a counterexample. Section [6] presents algorithms for nested sparsity using peb- 
ble games as oracles. Finally, Section [7] discusses extensions, applications and 
future directions. 

2 Preliminaries 

Our results rely on the same mathematical background as the work on body-and- 
bar rigidity by Tay [301 and White and Whiteley [34] . We use Grassmann-Cayley 
algebra, Plucker coordinates and instantaneous screw theory (see, e.g., [32] [33] 
and [22]). For self-containment, we briefly introduce notation and basic concepts 
from the Grassmann-Cayley algebra and its correspondence with instantaneous 
screws. 

2.1 Terminology and notation. 

We restrict ourselves to dimension 3 in this paper; 2-tensors in the Grassmann- 
Cayley algebra (see, e.g., [32] [33]) are identified with vectors in R 6 . The 
Grassmann-Cayley join operator is represented with V. The join p V q of two 
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vectors p, q G M 4 is the collection of all 6 minors of the matrix M obtained with 
p and q as its rows. We fix a convention at this point to order the minors in a 
6-vector as (|Mi 4 |, |M 24 |, |M 34 |, |Af 23 |, -|Mi 3 |, \M 12 \f\ The dot product of two 
vectors u and v is denoted (u, v). The star operator * swaps the first and last 3 
coordinates of a 6-vector. If p G M 3 and c G R, we denote by (p : c) the vector 
of length 4 obtained by appending c to p. 

Rigid body motions. The theory of screws was introduced by Ball pQ as a 

way of expressing rigid body motion. Rigid body transformations are associated 
with elements of the special Euclidean group SE(3). By Chasles' Theorem 
from 1830 (see [22]), they can also be expressed as screw motions (see Figure 



la). It follows that every instantaneous rigid body motion can be expressed 



as an instantaneous screw motion (see Figure [Tb|) ; for further details, we refer 
the reader to a standard text, e.g., page 24 of [22 . Both screw motions and 



(a) Every rigid body 
motion can be expressed 
as a screw motion that 
includes rotation and 
translation along the 
screw axis. 



(b) An instantaneous 
screw has both rota- 
tional and translational 
components. 



Figure 1: A screw motion and its associated instantaneous screw motion. 



instantaneous screw motions are defined with respect to a screw axis along with 
a rotation about the axis and a translation along it. 

In this paper, we are concerned only with instantaneous screw motions, 
which, for brevity, will be referred to as instantaneous screws. An instantaneous 
screw is represented by a 6-vector s = (— u?,v), where u?,v G M 3 ; the minus 
sign in front of uj is a convenient, technical convention. The first component u> 
encodes the angular velocity; as a vector, uj gives the direction of the screw axis, 

^■We remark that other papers (e.g., 30 and 34 ), use a different convention by fixing the 
order as (|Afi 2 |, |Mi 3 |, |M M |, |M 23 |, |M 24 |, |M 34 |) 
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and its magnitude encodes the angular speed. The translational velocity can be 
computed from uj and v, but we skip the details as they are not relevant for the 
rest of the paper. Note that the star operator applied to a screw s = (— o;,v) 
gives s* = (v, — ui). There is an exact correspondence between 2-tensors and 
instantaneous screws. This correspondence is the key to describing the rigidity 
matrix of the body-and-cad structures. 

2.2 Body-and-cad structures 

A body-and-cad structure in 3D is composed of n bodies interconnected by pair- 
wise constraints. Each body i is represented by a frame of reference, specified by 
a transformation matrix T{ from the special Euclidean group SE(3). Each body 
i additionally has a set of geometric elements (points, lines or planes) identified 
as attachments for the constraints. 

Representation of geometric elements. Each geometric element is rigidly 
affixed to a body i and is described with coordinates that are local with respect 
to the frame of reference for body i. For ease of analysis, we represent a plane 
in point-normal form as the pair (p,d), with p, d G M 3 , where p is a point on 
the plane and d is the normal to the plane. We represent a line in parametric 
form, given by the pair (p,d), with p, d G M 3 , where p is a point on the line 
and d is its direction. 

2.2.1 Cad graphs 

We now introduce the cad graph, our main combinatorial object for body-and- 
cad rigidity. 

To illustrate this concept, consider the following example, depicted in Figure 
[2] Let A and B be two dice rigidly stacked with the following constraints: (i) 
(plane-plane parallel) A : s Face 1 is parallel to £?'s Face 1, (ii) (plane-plane 
perpendicular) A's Face 2 is perpendicular to £?'s Face 3, (hi) (line-plane 
distance) The distance between A's Line 12 (intersection of Faces 1 and 2) and 
B's Face 1 is 1, and (iv) (point-point coincidence) A's Corner 236 (the point 
defined by Faces 2, 3 and 6) is coincident to £?'s Corner 123. These constraints 
are captured by a graph with two nodes connected by annotated edges called 
the cad graph; see Figure [3j 

Formally, a cad graph (G, c) is a multigraph G = ( V, E) together with an 
edge coloring function c : E — » C, where C = {ci, C2, . . . , C21} consists of 21 
colors corresponding to the full set of cad constraints: 

1. point-point coincidence 

2. point-point distance 

3. point-line coincidence 

4. point-line distance 

5. point-plane coincidence 
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Figure 2: Two dice rigidly stacked; die A is above B. Faces are labeled by the 
number of dots, and face 6 lies at the bottom (opposite 1). The length of an 
edge is 1. 




Figure 3: The cad graph for the example depicted in Figure 2 



6. 


point-plane distance 


7. 


line-line parallel 


8. 


line-line perpendicular 


9. 


line-line fixed angular 


10. 


line-line coincidence 


11. 


line-line distance 


12. 


line-plane parallel 


13. 


line-plane perpendicular 


14. 


line-plane fixed angular 


15. 


line-plane coincidence 



7 



16. 


line-plane distance 


17. 


plane-plane 


parallel 


18. 


plane-plane 


perpendicular 


19. 


plane-plane 


fixed angular 


20. 


plane-plane 


coincidence 


21. 


plane-plane 


distance 



The geometric meaning of these colors will be described in the next section; the 
cad graph only captures the type of constraint imposed. 

2.2.2 Body- and- cad frameworks 

A body-and-cad framework (G, c, L±, . . . , L21) is a cad graph (G, c) along with a 
family of functions Li, . . . , L21 describing the geometry of the structure, where 
the function Li captures the constraints corresponding to edges with the zth 
color. Let Ei = {e G E\c(e) = q} be the set of Q-colored edges. 

For example, the function for plane-plane fixed angular constraints L19 : 
£19 ->> (R 3 xR 3 )x(R 3 xR 3 )xR maps an edge e = zj to a triple ((p^, d^), (p^, dj), a) 
so that the planes (p^,d^) and (p^dj) affixed to bodies z, respectively j, are 
constrained to have the angle a between them. The function for point-point 
coincidence constraints L\ : E\ — » R 3 x R 3 maps an edge e = zj to a pair of 
points (pi,Pj) affixed to bodies z, respectively j, that are constrained to be co- 
incident. We will define the complete family of functions when analyzing them 
individually in Section [4] 

A realization G(T) of a body-and-cad framework (G, c, Li, . . . , L21) assigns 
a tuple of frames T = (Ti,...,T n ) for each vertex, satisfying the specified 
constraints. In this paper, we are not concerned with realization questions. 
We will always assume that a body-and-cad framework is given by a concrete 
realization, from which the family of functions Li, . . . , L21 is computed. 

Body-and-cad rigidity. Intuitively, a body-and-cad framework is rigid if the 
only motions respecting the constraints are the trivial 3D motions (rotations 
and translations); otherwise, it is flexible. We omit the technical definition, as 
it falls outside the scope of this paper. 

Body-and-cad minimal rigidity. For classical distance constraints, the con- 
cept of minimal rigidity is defined as follows: a structure is minimally rigid if the 
removal of any constraint results in a flexible structure. However, in our case, 
geometric constraints may correspond to more than one "primitive" constraint. 
Intuitively, a primitive constraint yields only one row in the rigidity matrix 
(formally defined in Section [3J, while the body-and-cad constraints may yield 
several rows. In our setting, we define minimal rigidity as above, but referring 
to the removal of primitive constraints only: a rigid body-and-cad structure is 
minimally rigid if the removal of any primitive constraint results in a flexible 
structure. 
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We return to the example from Figure [2] to illustrate the subtleties of this 
concept. The structure depicted is rigid. We say the structure is over constrained 
since it remains rigid even after the removal of constraint (iii). The resulting 
structure is now minimally rigid. As we will see in Section [3j constraints (i), 
(ii) and (iv) correspond to 6 primitive constraints. Thus, the removal of any 
primitive constraint results in a flexible structure. 

Now consider stacking the dice with the following two constraints: (i) (line- 
line coincidence) A's Line 26 is coincident to £?'s Line 12 and (ii) (line-line 
coincidence) A : s Line 36 is coincident to £?'s Line 13. This structure is still 
rigid. While it becomes flexible after the removal of either constraint (i) or (ii), 
it is not minimally rigid. As we will see in Section [3J a line-line coincidence 
constraint corresponds to 4 primitive constraints. Thus, this structure has 8 
primitive constraints and is over constrained. To give some intuition, note that 
a structure composed of 2 rigid bodies has 12 degrees of freedom. Of these, 6 
are trivial, so we may fix body A to factor them out. Now consider constraint 
(i); the structure is left with 2 degrees of freedom, as B may slide along the 
line and rotate about it. This line-line coincidence constraint is "eliminating" 
4 degrees of freedom, formalized by the 4 rows of the rigidity matrix developed 
in Section [3] for the line-line coincidence constraint. 

Body-and-cad infinitesimal rigidity. Infinitesimal rigidity is the linearized 
version of rigidity and is the only type we study in this paper. Let s = 
(si,...,s n ) G (M 6 ) n assign an instantaneous screw to each body i and let 
s* = (s*,...,s*). The vector s is an infinitesimal motion of a body-and-cad 
structure if it infinitesimally respects the constraints. This can be expressed 
with the help of the rigidity matrix, fully described in Section [3j An infinitesi- 
mal motion is a vector in the kernel of the rigidity matrix. The kernel always 
contains the trivial infinitesimal motions, defined as those s with = Sj for all 
i and j. 

A body-and-cad framework is infinitesimally rigid if the only infinitesimal 
motions are trivial; otherwise, it is infinitesimally flexible. 

Remarks. To develop the rigidity theory for a new model, three steps must be 
accomplished. 

1. Algebraic theory. Formulate the rigidity concept in algebraic terms, re- 
sulting in an algebraic variety. 

2. Infinitesimal theory. Analyze the local behavior at some point on the 
algebraic variety. This reduces to the study of a rigidity matrix. 

3. Combinatorial rigidity. Seek a combinatorial characterization of minimal 
rigidity in terms of properties of an underlying graph structure. This is 
usually derived from properties of the rigidity matrix at a generic point 
on the algebraic variety. 

In this paper, we directly formulate the infinitesimal rigidity theory for body- 
and-cad structures and identify combinatorial properties for the generic case. 



9 



To summarize, a cad graph is an edge- colored multigraph that captures the 
body-and-cad combinatorics, and a body- and- cad framework captures the geom- 
etry of the structure. As we develop the analysis of these concepts, an additional 
combinatorial object called the primitive cad graph will be associated to the cad 
graph. This is a multigraph with red or black edges, which captures certain 
combinatorial properties of infinitesimal body-and-cad rigidity. 



3 Foundations of infinitesimal theory 

The example from the previous section exposes some of the subtleties encoun- 
tered with body-and-cad constraints that are not found when considering classi- 
cal distance constraints. We introduce two new concepts to simplify the analysis: 
primitive angular and blind constraints. We then define, as building blocks, 4 
basic angular and blind constraints and develop their infinitesimal theory. All 21 
body-and-cad constraints can be studied using these building blocks, leading 
to the body-and-cad rigidity matrix. 



3.1 Primitive constraints 

A primitive constraint is one that may restrict at most one degree of freedom. 
For example, a point-point distance (bar) constraint is a primitive constraint, 
while the line-line coincidence constraint from the example in the preceding 
section is not. We classify primitive constraints into two types: angular and 
blind; as the theory is developed, it will become more clear why these classifica- 
tions are appropriate, as they correspond to constraints demonstrating different 
algebraic behaviors. 

A rigid body in 3D has 6 degrees of freedom, 3 of which are rotational and 
3 of which are translational. A primitive angular constraint may restrict only a 
rotational degree of freedom, whereas a primitive blind constraint may restrict 
either a rotational or a translational degree of freedom. For instance, a line-line 
perpendicular constraint is a primitive angular constraint as it may restrict 
at most one rotational degree of freedom. A point-point distance (bar) con- 
straint is a primitive blind constraint as it may restrict at most one rotational 
or translational degree of freedom. We will associate a set of primitive angular 
and a set of primitive blind constraints with each body-and-cad constraint. 



3.2 Rigidity matrix 

The rigidity matrix R for a body-and-cad structure has 6 columns for each body 
z, corresponding to the components of the instantaneous screw s^, as was done 
for the original body-and-bar rigidity matrhQ There is a row for each primi- 
tive constraint associated to the original body-and-cad structure. A primitive 
angular constraint results in a row containing zero entries in the first 3 columns 



2 The starred version s* (see Section 2.1 ) will be used to conveniently order the columns of 
the rigidity matrix. 
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for each body, while a primitive blind constraint may have non-zero entries in 
any of the 6 columns for each body. In the schematic below, gray cells indicate 
potentially non-zero entries, and red cells highlight the zero entries for angular 
constraints. 

Since the trivial motions corresponding to the 3D rigid motions are neces- 
sarily in the kernel of R, the maximum rank of R is 6n — 6. By definition, a 
structure is infinite simally rigid if its rigidity matrix has rank exactly 6n — 6. 
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3.3 Building blocks 

We now define 4 very specific basic angular and blind constraints (2 of each) and 
develop the infinitesimal theory for them. Everything in Section [4] is derived 
from these basic building blocks. The material presented here is the most 
technical part of our paper. 
Angular building blocks 

All body-and-cad angular constraints can be reduced to the following basic 
constraints between pairs of lines: 

(i) basic line-line non-parallel fixed angular, and 

(ii) basic line-line parallel. 

(i) Basic line-line non-parallel fixed angular. A line-line non-parallel angular 
constraint between bodies i and j is defined by identifying a pair of non-parallel 
lines, each rigidly affixed to one body, and fixing the angle between them. Let 
di and dj be the directions of the lines affixed to bodies i and j, respectively. 
Then the constraint is infinitesimally maintained if the axis of the relative screw 
— Sj is in a direction lying in the plane determined by d^ and dj, i.e., 

{{uji -<A>j),di x dj) = 

Since — is composed of the last three coordinates of s*, this is equivalent to 

<(s*-s*),((0,0,0),d,- xd,)> = (3.1) 

This corresponds to one row in the rigidity matrix: 
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dj x di 


-0- 







-0- 



i?asic line-line parallel constraint. A line-line parallel constraint between 
bodies z and j is defined by identifying a pair of parallel lines, each rigidly 
affixed to one body, and restricting them to remain parallel. Let d = (a, 6, c) 
be the direction of the parallel lines. Then the constraint is infinitesimally 
maintained if the axis of the relative screw — Sj is in the same direction as d, 
i.e., (u>i — ojj) = ad, for some scalar a. This can be expressed by the following 
two linear equations, where uj = u)i — ujj = (uo x ,uj y ,uj z ): 



uj y c-uj z b 



Since — is composed of the last three coordinates of s*, these are equivalent 
to 



((s*-s*),(0,0,0,-M,0)) 
((s*-s*),(0,0,0,0,-c,6)) 

and correspond to two rows in the rigidity matrix: 



(3.2) 
(3.3) 



•0- 



(-M,o) 


...0- 




(6,-0,0) 


...o... 


(0,-c,b) 


-0- 




(0,c,-b) 


...Q.. 



Blind building blocks 

Let p be a point and p' its instantaneous velocity resulting from the instanta- 
neous screw s G M 6 . Let c G I 3 be an arbitrary direction vector. We either 
constrain the velocity p' to be orthogonal or parallel to c. This yields the re- 
maining basic constraints: 



(Hi) basic blind orthogonality (see Figure 4a), and 
(iv) basic blind parallel (see Figure 4b). 

Expressing both of them becomes straightforward using the following fact: 



Fact 1. Let s G M 6 be an instantaneous screw, p G 
velocity of p under the screw motion s. Then 

sV(p:l) = (p / ,-(p,p / )) 

and, for any q G R 3 and q w G R, 

sV(p:l)V(q:0 = (p / ,q>-f <P,p'> 



a point and p' the 
(3.4) 

(3.5) 



Proof In the following, superscripts x,y,z,w denote the components of a vector 
in R 4 . The minor of a 3 x 4 matrix A determined by columns z, j and k is denoted 
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If s is a decomposable 2-tensor (a 2-extensor), then its components are the 
minors of a 2 x 4 matrix M; see, e.g., [32] [33] for a standard review of 2-tensors 
in Grassmann-Cayley algebra. Let A be the 3x4 matrix obtained by appending 
(p : 1) to the bottom of M. The join sV(p : 1) is the collection of the four minors 



of A. We fix the convention that s V (p : 1) 
Then 



s V (p : 1) V (q : q v 



^234 1 



M 



1134| 



24 | 



|Ai 23 |). 



q z 



Performing a Laplace expansion along the 4th row of the matrix yields q x | A234I — 
q y \AM + 9l^i2 4 | - q w \A 12 s\ = q x (s V (p : l)) x + ^(s V (p : 1))* + q z (s V (p : 
l)) z +f(sV(p: l)) w . 

Crapo and Whiteley [2 derived that s V (p : 1) = (p' : — (p, p')). Applying 
it, we obtain our desired result. The derivation when s is indecomposable (the 
sum of two 2-extensors) is a simple extension obtained by working with the two 
2-extensors simultaneously. □ 





(a) Orthogonality: constrains the instanta- (b) Parallel: constrains the instantaneous 

neous velocity p ; of the point p to be or- velocity p ; of a point p to lie in the same 

thogonal to a direction c. Then p x must lie direction as a vector c. 
in the plane (p, c). 

Figure 4: Basic blind geometric constraints. 



(Hi) Basic blind orthogonality constraint 

This constrains the instantaneous velocity p' of the point p to be orthogonal to 
a direction c. To express this (see Figure 4a), we simply substitute q = c and 



into Equation 3.5 Then (p',c) = if and only if 
s V (p : 1) V (c : 0) = 



if and only if 



(s*,(p:l)V(c:0))=0 



(3.6) 



(iv) Basic blind parallel constraint 

This constrains the instantaneous velocity p' of the point p to lie in the same 



direction as a direction c. To express this (see Figure 4b), we apply Equation 
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3.5 twice by substituting q = (c y , — c x , 0) and q w = first, then q = (0, c 2 , - 



and <2 W = 0. We obtain that p' = ac for some a G R if and only if 



if and only if 



sV(p: 1) V(c^,-c*,0,0) = 
sV(p:l)V(0,c^,-c y ,0) = 



<s*,(p:l)V(^,-c*,0,0)) = 
<s*,(p:l)V(0,c*,-^,0)> = 



(3.7) 
(3.8) 



4 Infinitesimal theory for body- and- cad constraints 

We use the four basic building blocks just presented to complete the development 
of the infinitesimal theory. In this section, we present the rows of the rigidity 
matrix associated with each of the 21 body-and-cad constraints. In all figures, 
body i is represented by the green tetrahedron and body j by the purple cube. 



4.1 Angular constraints 

Angular constraints may be parallel, perpendicular or arbitrary fixed angular 
constraints; see Figures [5}|7| for depictions of line-line, line-plane, and plane- 
plane angular constraints. 






(a) Parallel. (b) Perpendicular. (c) Arbitrary fixed angular. 

Figure 5: Line-line angular constraints. 

We now formally define the functions for describing angular constraints. 

• Line-line parallel: Lj : Ej — >> R 3 x R 3 x R 3 maps an edge e = ij to 
a triple (pi,Pj,d) so that the lines (pi,d) and (pj,d) affixed to bodies i 
and j, respectively, are constrained to remain parallel to each other. 
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(a) Parallel. (b) Perpendicular. (c) Arbitrary fixed angular. 



Figure 6: Line-plane angular constraints. 




(a) Parallel. (b) Perpendicular. (c) Arbitrary fixed angular. 



Figure 7: Plane-plane angular constraints. 

• Line-line perpendicular: L 8 : Eg —> (R 3 x R 3 ) x (R 3 x R 3 ) maps an 
edge e = ij to a pair ((pi, d^), (pj, dj)) so that the lines (p^,d^) and 
(pj,dj) affixed to bodies i and j, respectively, are constrained to remain 
perpendicular to each other. 

• Line-line fixed angular: L 9 : E 9 (R 3 x R 3 ) x (R 3 x R 3 ) x R maps an 
edge e = ij to a triple ((p^d^), (pj,dj),a) so that the lines (p^d^) and 
(pj, dj) affixed to bodies i and j, respectively, are constrained to maintain 
the angle a between them. 

• Line-plane parallel: L 12 : E 12 — > (R 3 x R 3 ) x (R 3 x R 3 ) maps an edge 
e = ij to a pair ((pi, d^), (pj, dj)) so that the line (p^,d^) and plane 
(pj,dj) affixed to bodies i and j, respectively, are constrained to remain 
parallel to each other. 

• Line-plane perpendicular: Lis ' E13 -> (^ 3 x ^ 3 ) x (^ 3 x ^ 3 ) maps an 
edge e = ij to a pair ((p^, d^), (pj, dj)) so that the line (p^, d^) and plane 
(pj,dj) affixed to bodies i and j, respectively, are constrained to remain 
perpendicular to each other. 

• Line-plane fixed angular: L 14 : E 14 (R 3 x R 3 ) x (R 3 x R 3 ) x R maps 
an edge e — ij to a triple ((pi, d$), (pj, dj), a) so that the line (p^, d^) and 
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plane (pj,dj) affixed to bodies i and j, respectively, are constrained to 
maintain the angle a between them. 

• Plane-plane parallel: L17 : £17 — ^ R 3 x R 3 x R 3 maps an edge e = ij to 
a triple (p^,pj,d) so that the planes (pi,d) and (pj,d) affixed to bodies 
i and j, respectively, are constrained to remain parallel to each other. 

• Plane-plane perpendicular: L 18 : Eis — » (R 3 x R 3 ) x (R 3 x R 3 ) maps 
an edge e = ij to a pair ((pi, d^), (pj, dj)) so that the planes (p^, d^) and 
(pj,dj) affixed to bodies i and j, respectively, are constrained to remain 
perpendicular to each other. 

• Plane-plane fixed angular: L 19 : £19 (R 3 x R 3 ) x (R 3 x R 3 ) x R 

maps an edge e = ij to a triple ((pi, d^), (p^, dj), a) so that the planes 
(p^ d^) and (pj, dj) affixed to bodies i and j, respectively, are constrained 
to maintain the angle a between them. 

It is straightforward that the line-line angular constraints (perpendicular, 
fixed angular and parallel) are expressed by the two basic angular building 
blocks. We observe that the line-plane and plane-plane angular constraints 
reduce to them as follows. 

• Line-plane parallel: Reduces to line-line non-parallel fixed angular 

using the normal to the plane. 

• Line-plane perpendicular: Reduces to line-line parallel using the 
normal to the plane. 

• Line-plane fixed angular: Reduces to line-line non-parallel fixed 
angular using the normal to the plane. 

• Plane-plane parallel: Reduces to line-line parallel using the planes' 
normal. 

• Plane-plane perpendicular: Reduces to line-line non-parallel fixed 
angular using the planes' normals. 

• Plane-plane fixed angular: Reduces to line-line non-parallel fixed 
angular using the planes' normals. 



4.2 Blind constraints 

The remaining coincidence and distance constraints reduce to some combination 
of basic angular and basic blind constraints from Section 3.3 We consider them 
in the following order: point-point, point-line, point-plane, line-line, line- 
plane and plane-plane. Since a point-point distance constraint (Figure 
[#&[ )) is a bar (see [30l [34]), we consider only the point-point coincidence 
constraint. 

The function L\ : E\ — >> R 3 maps an edge e = ij to a point p so that it is 
constrained to lie on bodies i and j simultaneously. 
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(a) Coincidence: bodies i and j must coin- 
cide on the specified point p. 



(b) Distance (bar): a point affixed to body 
i must be a fixed distance from a point af- 
fixed to body j. 



Figure 8: Point-point constraints. 

Then the point-point coincidence constraint (Figure\8afy is infinitesimally 
maintained if the relative velocity of p is exactly 0. Since = (— u^, v^) and Sj = 
(— ojj, Vj), then the relative screw is defined by (— u>, v), where uj = {ui — u)j) 
and v = (v^ — Vj). Therefore, the constraint is infinitesimally maintained if and 
only if p's infinitesimal velocity p' = wxp + v = 0, i.e., 

wY-wy + u* = o 

uj z p x - uj x p z + v y = 
uj x p y - uj y p x + v z = 



if and only if 



<(*; 

{{s*- s*), (0,0,1,-^,^,0)) = 



s*), (1,0, 0,0, -p z ,p y )) 
s*), (0,1, 0,p z ,0,-p x )) 



(4.1) 
(4.2) 
(4.3) 



Thus, a point-point coincidence constraint corresponds to 3 rows in the rigidity 
matrix: 



a. 





(1,0,0,0,-^,^) 




(-1,0,0,0,^,-^) 


...o... 


-o- 


(0,1,0,^,0,-^) 


...o... 


(0,-1,0,-^,0,^) 


...O- 


-0- 


(0,0,1,-^,^,0) 


-0- 


(0,0,-1,^,-^,0) 





Point-line coincidence (Figure 9a): The function L 3 : E 3 — >> R 3 x (R 3 x R 3 ) 
maps an edge e = zj to a pair (p^, (pj, d)) so that point p^ affixed to body i is 
constrained to lie on the line (pj, d) affixed to body j. 

The point-line coincidence constraint is infinitesimally maintained by us- 
ing 2 primitive blind constraints from Equations |3.7| and |3.8| to express that the 
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(a) Coincidence: a point affixed to body (b) Distance: a point affixed to body i 
i must lie on a line with direction d affixed must be a fixed distance from a line affixed 
to body j. Then the instantaneous velocity to body j. If d is the perpendicular direc- 
p^ of pi must lie in the same direction as tion from the line to p^, then the instanta- 
d. neous velocity p^ of p^ must be orthogonal 

to d. I.e., p^ must lie in the plane (pj,d). 

Figure 9: Point-line constraints. 



relative velocity of lies in the same direction as d: 

((s i -s J )*,( Pi :l)v(^,-^,0,0)) 
((s i -s i )*,(p i :l)V(0,#,-^,0)) 



(4.4) 
(4.5) 



Thus, a point-line coincidence constraint corresponds to 2 rows in the rigidity 
matrix: 



-Mi 



...O- 


(Pi : 1) V(<jy,-(F, 0,0) 


...O- 


-(( Pi :l)V(^,-d a ,0,0)) 


...O- 


-0- 


( Pi : 1) v (0,^,-^,0) 


-0- 


-((Piii)v (0,^,-^,0)) 


-0- 



Point-line distance (Figure\9b\): The function L 4 : E A W x (R 3 x M 6 ) x K 
maps an edge e — ij to a triple (p^, (pj, d), a) so that point p^ affixed to body 
i is constrained to lie a distance a from the line (pj, d) affixed to body j. 

Let d be the perpendicular direction from the line (pj,d) to p^. Then the 
point-line distance constraint is infmitesimaliy maintained using 1 primitive 
blind constraint from Equation |3.6| to express that the relative velocity of p^ is 
orthogonal to d: 

'(s,-s,)*,(Pz : 1) V(d:0)) =0 (4.6) 



Thus, a point-line distance constraint corresponds to one row in the rigidity 
matrix: 
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-0- 


(Pi : 1) V (d : 0) 


•0- 


-(( Pi : 1) V (d : 0)) 


•••0- 





id) 


< 


> 


d 













(a) Coincidence: a point p-i affixed to body 
i must lie in a plane with normal d affixed 
to body j. Then the instantaneous velocity 
of pi must remain in the plane (p^,d). 



(b) Distance: a point p^ affixed to body i 
must be a fixed distance from a plane with 
normal d affixed to body j. Then the in- 
stantaneous velocity p^ of pi must remain 
in the plane (p^,d). 



Figure 10: Point-plane constraints. 



Point-plane coincidence (Figure 10a): The function L5 : E$ 



maps an edge e = ij to a pair (p^, (pj, d)) so that the point p^ affixed to body 
i is constrained to lie in the plane (pj, d) affixed to body j. 

The point-plane coincidence constraint is infinitesimally maintained us- 



ing Equation 3.6 to express that the relative velocity of p^ remains in the plane: 

((8 i -8 j )%(Pi:l)V(d:0)>=0 (4.7) 

Thus, a point-plane coincidence constraint corresponds to one row in the rigidity 
matrix: 



0- ( Pi : 1) V (d : 0) 



•0" 



-(( Pi : 1) V (d : 0)) 



Point-plane distance (Figure 10b): The function L 6 : E 6 — > R 3 x (R 3 

maps an edge e = ij to a triple (p^, ( Pj -,d),a) so that the point p^ affixed to 
body i is constrained to lie a distance a from the plane (pj , d) affixed to body 
j- 

The point-plane distance constraint is infinitesimally maintained by using 
Equation |3.6| to express that the relative velocity of p^ remains parallel to the 
plane: 



<(s,- Sj )*,(Pz : 1) V(d:0)> =0 



(4.S 
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Thus, a point-plane coincidence constraint corresponds to one row in the rigidity 
matrix: 



0- ( Pi : 1) V (d : 0) -0- -(( Pi : 1) V (d : 0)) -0 




(a) Coincidence: bodies i and j must coin- 
cide on the specified line (p, d). Then, in 
addition to a line-line parallel constraint, 
the instantaneous velocity p ; must lie in 
the same direction as d. 




(b) Distance: a line (p^, d^) affixed to body 
i must be a fixed distance from a line 
(pj, dj) affixed to body j. If p is the point 
on the line (p^, d^) that is closest to the line 
(pj, dj), then the instantaneous velocity p x 
of p must be orthogonal to d^ x dj , the per- 
pendicular to both lines. I.e., p' must lie 
in the plane (p,d^ x dj) 



Figure 11: Line-line constraints. 



3 x M 3 maps 



Line-line coincidence (Figure 11a): The function Lio : Eio — » 
an edge e = ij to a pair (p, d) so that the line (p, d) is constrained to be affixed 
to bodies i and j simultaneously. 

We place a line- line parallel angular constraint, resulting in 2 primitive 
angular constraints from Equations |3.2| and |3.3| 



<(s*-s*), (0,0, 0,-dV*,0)) (4.9) 
((s*-s*), (0,0,0,0,-^)) (4.10) 

Then, to maintain coincidence, associate 2 primitive blind constraints from 



Equations 3.7 and|3.8|to force the relative velocity of p to lie along d: 



(( Si - Sj )*,(p:l)V(^,-^,0,0)) = 
(( Si - Sj )*,(p:l)V (0,<f,-^,0)) = 



(4.11) 
(4.12) 
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These 4 equations maintain the line-line coincidence constraint infinites- 
imally and correspond to 4 rows in the rigidity matrix: 















...o... 





(-dy,d x ,0) 




"on (^,-^,o) 


...o... 


...o... 




(0,-d z ,dy) 




~o\ (o,d z ,-dy) 


...o... 


...o... 


(p: 


1) V (^,-^,0,0) 


...o... 


-((p:l)V(d^-d*,0,0)) 


...o... 


-0- 


(P: 


1) V (0,^,-^,0) 


-0- 


-((p:l)V (0,^,-^,0)) 


-0- 



Line-line distance (Figure ] lib]) : The function L n : En (R 3 x R 3 ) x (R 3 x 
R 3 ) x R maps an edge e = zj to a triple ((p^d^), (pj,dj),a) so that the lines 
(pi,di) and (pj,dj) affixed to bodies i and j, respectively, are constrained to 
lie a distance a from each other. 

Let p G R 3 be the point on the line (p^, d$) closest to the line (pj, dj). Then 
the line-line distance constraint is infmitesimaHy maintained if the relative 
velocity of p is orthogonal to the direction perpendicular to both lines. In other 
words, p' must lie in the plane defined by the point p and normal direction 
di x dj. By substituting p and d^ x dj into Equation |3.6[ we obtain the linear 
equation 



((s i -s J )*,( P :l)V((d i xd i ):0)) = 



(4.13) 



associating one blind primitive constraint. This corresponds to one row in the 
rigidity matrix: 



(p : 1) V ((di x d 3 ) : 0) -0- -((p : 1) V ((d< x d,-) : 0)) 





(a) Coincidence: a line (pi,dj) affixed to 
body i must lie in a plane with normal dj 
affixed to body j. Then, in addition to a 
line-plane parallel constraint, the instanta- 
neous velocity of must lie in the plane 
(Pijdj)- 



(b) Distance: a line (p^, d^) affixed to body 
i must be a fixed distance from a plane with 
normal dj affixed to body j. Then, in addi- 
tion to a line-plane parallel constraint, the 
instantaneous velocity p^ of pi must lie in 
the plane (p^, dj). 



Figure 12: Line-plane constraints. 
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Line-plane coincidence (Figure 12a): The function L15 : £15 — (R x R 3 ) x 

(R 3 x R 3 ) maps an edge e = ij to a pair ((p^, d;), (pj, dj)) so that the line 
(pi,di) affixed to body i is constrained to lie in the plane (pj,dj) affixed to 
body j. 

The line-plane coincidence constraint is infmitesimaHy maintained using 



a primitive angular line-plane parallel constraint from Equation 3.1 
((s*-s*), ((0,0,0),^ xd,)) = 



(4.14) 



In addition, a primitive blind constraint from Equation 3.6 forces the relative 
velocity of p^ to remain in the plane: 



((s i -s J )*,( Pi :l)V(d i :0))=0 



(4.15) 



Thus, a line-plane coincidence constraint corresponds to 2 rows in the rigidity 
matrix: 



•0" 



dj x di 



-0- ( Pi : 1) V (dj : 0) 



•0- 



d ? x d, 



-0- -(( P< : 1) V (d, : 0)) 



Line-plane distance (Figure 12b): The function Liq : Eiq 



(R 3 x R 3 ) x R maps an edge e = ij to a triple (( P ;, d^), (pj, dj), a) so that the 
line (p^d^) affixed to body i is constrained to lie a distance a from the plane 
(pj,dj) affixed to body j. 

The line-plane distance constraint is maintained infinitesimally by using 
the same equations as for the line-plane coincidence constraint: a primitive 



angular line-plane parallel constraint from Equation 3.1 



((s*-s*), ((0,0,0),^ xd,)) = 



(4.16) 



In addition, a primitive blind constraint from Equation 3.6 forces the relative 
velocity of p^ to remain parallel to the plane: 



<(s i -s J -r,(Pi:l)V(d j :0))=0 



(4.17) 



Thus, a line-plane distance constraint corresponds to 2 rows in the rigidity 
matrix: 



•0" 



dj x di 



-0- ( P< : 1) V (dj : 0) 



d ? x d, 



•0" 



-((Pi : 1) V (d, : 0)) 



•0- 



Plane-plane coincidence (Figure 13a): The function L20 : £^20 ^ R x R 3 
maps an edge e = ij to the pair (p, d) so that the plane (p, d) is constrained to 
be affixed to both bodies i and j simultaneously. 
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(a) Coincidence: bodies i and j must co- 
incide on the specified plane (p, d). Then, 
in addition to a plane-plane parallel con- 
straint, the instantaneous velocity p ; of p 
must remain in the plane. 



(b) Distance: a plane (p^,d) affixed to 
body i must be a fixed distance from a 
plane (pj,d) affixed to body j. Then, 
in addition to a plane-plane parallel con- 
straint, the instantaneous velocity p^ of p^ 
must remain in the plane (p^,d). 



Figure 13: Plane-plane constraints. 



We place a plane-plane parallel angular constraint, resulting in 2 primitive 



angular constraints from Equations |3.2| and 3.3 



((s*-s*), (0,0, 0,-d?,<P,0)) 
<(s*-s*),(0,0,0,0,-eF,^)) 



(4.18) 
(4.19) 



Then, to maintain coincidence, place a primitive blind constraint by using Equa- 
tion [3]6] to force the relative velocity of p to remain in the plane: 

<(s,- Sj )*,(p: 1) V(d:0)> =0 (4.20) 
Thus, a plane-plane coincidence corresponds to 3 rows in the rigidity matrix: 



-0- 1 1 (-d y ,d x ,0) 


...o- 




((P,-d x ,0) 




-0- 1 


(p,-d z ,dy) 


...o... 





(0,d*,-d*) 


...o... 


-0- (p : 1) V (d : 0) 


-0- 


-((p : 1) V (d : 0)) 


-0- 



Plane-plane distance (Figure] 13b]/: The function L 2i : E 2 i M 3 xR 3 xR s xl 
maps an edge e = ij to a quadruple (p^,pj,d,a) so that the planes (pi,d) 
and (pj,d) affixed to bodies i and j, respectively, are constrained to have the 
distance a between them. 

Maintaining the plane-plane distance innnitesimally reduces to the same lin- 
ear equations as for plane-plane coincidence. We place a plane-plane parallel 
angular constraint, resulting in 2 primitive angular constraints from Equations 



3.2 and 3.3 along with a primitive blind constraint using Equation 3.6 to force 



23 



the relative velocity of to remain parallel to the plane: 
((s*-s*), (0,0, 0,-^,^,0)) 



)), (0,0,0,0, -d',^)) 



(( Si - Sj )*,( Pi : 1) v(d:0)> =0 



(4.21) 
(4.22) 

(4.23) 



Thus, a plane-plane distance constraint corresponds to 3 rows in the rigidity 
matrix: 



< 








Vj -U>i 








-0- (-d^d^O) 




0| (d»,-d*,0) 


...o... 


-0- | (0,-d z ,d y ) 




~cn (o.d'.-d") 


...O- 


-0- ( Pi : 1) V (d : 0) 


-0- 


-(( Pi : 1) V (d : 0)) 


-0- 



4.2.1 Example 

To help the reader, we complete the formalization of the dice example depicted 
in Figures [2] and [3] We assume that the z-axis lies along the base of the dice in 
the direction of Face 2, with the x?/-plane parallel to Face 3. 
Then the framework is described by the functions: 

• i 17 (e (i) ) = ((0,2,0), (0,1,0), (0,1,0)) 

• Li 8 (e (ii) ) = (((0, 2, 0), (1, 0, 0)), ((0, 0, 1), (0, 0, 1))) 

• L 16 (e (iii) ) = (((0, 2, 0), (0, 0, 1)), ((0, 1, 0), (0, 1, 0)), 1) 

• £i(e(»„)) = (0,1,1) 

Since the example only uses four types of constraints, we omit the description 
of the remaining Lj functions. 

For each edge, we construct the associated rows in the rigidity matrix, re- 
sulting in the following: 



e ( «) 











-LO A 
















-1 

































I 













1 


; 


1 


1 








o 


J 
















-1 

















1 





1 











-1 


1 


-1 











1 





1 











-1 











1 


-1 














-1 



-UJ B 


1 














-1 





-1 





-1 




















1 


-1 


-1 








1 
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4.3 Summary of infinitesimal theory 

We have now completed the development of the infinitesimal theory for body- 
and-cad rigidity. Table [I] summarizes the associations for each constraint to the 
number of primitive angular and blind constraints. As an example of how to 
read the table, the last two columns (corresponding to plane) of row 3 (corre- 
sponding to coincidence under line) indicate that a line-plane coincidence 
constraint reduces to 1 angular and 1 blind primitive constraint. In the next 
section, we identify a combinatorial property based on the shape of the rigidity 
matrix. 





point 
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angular 


blind 


angular 
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angular 
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distance 
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fixed angular 






1 
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2 


1 
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1 


parallel 










2 





perpendicular 










1 





fixed angular 










1 






Table 1: Association of body-and-cad (coincidence, angular, distance) con- 
straints with the number of blind and angular primitive constraints. 



5 Combinatorics 

Now we address the question of combinatorially characterizing when a body- 
and-cad rigidity matrix is generically independent, i.e., the rank function drops 
only on a measure-zero set of possible entries. The shape of the rigidity ma- 
trix leads to a natural property that we call nested sparsity. We show that 
nested sparsity is a necessary condition for body-and-cad rigidity and prove by 
a counterexample that it is insufficient. 

Nested sparsity. A graph on n vertices is (k,£) -sparse if every subset of n' 
vertices spans at most kn' — i edges; it is tight if, in addition, it spans kn — t 
total edges. 
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Let G = (V, RU B) be a graph with its edge set colored into red and black 
edges, corresponding to R and £?, respectively. We say that G is ^2,^2)- 
nested sparse if it is (fei, ^i) -sparse and Gr = (V, R) is (&2, -^-sparse; the graph 
is ki^'i) -nested tight if, in addition, G is (fei, ^i)-tight. Note that nested 

sparsity only makes sense when (&2, ^2)-sparsity is more restrictive than (fci, ^1)- 
sparsity. 

Primitive cad graphs. Given a cad graph (G, c), we define the primitive cad 
graph H = (V,RUB) to be the multigraph obtained by assigning vertices to 
bodies and constraints to disjoint edge sets R and 5, corresponding respectively 
to primitive angular and blind constraints. For each edge e with type q, asso- 
ciate primitive angular constraints to edges in R and primitive blind constraints 
to edges in B as described in Table [I] Figure [l4| depicts the primitive cad graph 
associated with the dice example from Figure [2j whose cad graph is depicted in 
Figure [3] 




Figure 14: The primitive cad graph for the example depicted in Figures [2] and 
[3j dashed edges denote R edge set. 

Theorem 5.1. Let H = (V, R U B) be the primitive cad graph associated to a 
body- and- cad framework, where R and B correspond to primitive angular and 
blind constraints, respectively. Then (6, 6, 3, 3) -nested tightness is a necessary 
condition for generic minimal body- and- cad rigidity. 

Proof. Let A be the rigidity matrix associated with G. Reorder the columns so 
that the first 3n columns correspond to the — uj elements of the screws and the 
last 3n columns correspond to the v elements. Reorder the rows to have the \R\ 
rows corresponding to primitive angular constraints first; since these rows have 
all Os in the last 3n columns, we simply consider the submatrix Ar defined by 
these \R\ rows and the first 3n columns. Then it is clear that (3, 0)-sparsity 
is necessary on Gr = (V, R). To see that (3, 3)-sparsity is necessary, we note 
that the 3-dimensional space of trivial motions of so(3) (infinitesimal rotations) 
is a subspace of the kernel of Ar. These are defined by the basis {pi,P2?P3}j 
where pi is the vector obtained by n copies of (1, 0, 0), P2 is the vector obtained 
by n copies of (0,1,0), and ps is the vector obtained by n copies of (0,0,1). 
Similarly, for the overall (6, 6)-sparsity, note that we have a 6-dimensional space 
of trivial motions of se(3) (infinitesimal rigid body motions), defined by the basis 
{Pi, P2, P3, Ti, T2, T3}, where pi simply appends 3n zeros to pi\ n is the vector 
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(a) Structure with 3 bodies and 6 constraints. A and B have 2 point-point distance 
constraints (denoted by the gray pair of points and black pair of points) and a line-line 
coincidence constraint (denoted by the shared dotted line) . A and C have a line-line fixed 
angular constraint (denoted by the dashed lines) and a plane-plane coincidence constraint 
(denoted by the two faces with vertical stripes). B and C have a line-plane coincidence 
constraint (denoted by the solid line on B and the checkered face on C). 





(b) The structure is flexible: C can move 
relative to A and B by translating in the 
direction indicated by the arrows. 



(c) Corresponding primitive cad graph is 
(6, 6, 3, 3)-nested tight; dashed edges de- 
note R edge set. 



Figure 15: Counterexample shows nested sparsity condition is not sufficient. 



obtained by n copies of (0, 0, 0, 1, 0, 0), T2 is the vector obtained by n copies of 
(0, 0, 0, 0, 1, 0), and T3 is the vector obtained by n copies of (0, 0, 0, 0, 0, 1). □ 

Counterexample. We now show that (6, 6, 3, 3)-nested sparsity is not sufficient 
for body-and-cad rigidity. The example in Figure [T5| depicts a flexible structure 
whose associated graph is (6, 6, 3, 3)-nested tight. It is composed of 3 bodies 
A, B and C; Figure [l5a| depicts the constraints. The structure has one degree of 
freedom, as indicated by the arrows in Figure 15b The associated primitive cad 



graph is shown in Figure 15c the reader may check that it is (6, 6, 3, 3)-nested 
tight. 



27 



6 Algorithms for nested sparsity 



In the previous section, we defined nested sparsity, proving that is a necessary, 
but not sufficient, condition for body-and-cad rigidity. We now examine the 
algorithmic aspects of nested sparsity. 




(a) Original graph is not (b) Subgraph that is max- (c) Subgraph that is max- 

(2, 2, 1, l)-nested sparse; imally (2, 2, 1, l)-nested imally (2, 2, 1, l)-nested 

dashed edges denote R edge sparse with 7 edges. sparse with 8 edges, 
set. 



Figure 16: Example showing why (2, 2, 1, l)-sparsity is not matroidal; two max- 
imal subgraphs do not have the same size. 



We first observe that nested sparsity is not matroidal, as seen by the exam- 



ple for (2, 2, 1, l)-nested sparsity in Figure 16 However, for certain values of 



&2 an d £2-, nested sparsity is the intersection of two matroids. 

Theorem 6.1. \3T$ When < £{ < 2ki, for i = 1,2, nested sparsity is the 
intersection of two matroids. 

Proof Define the ground set E to be the complete graph ^ 2/ci_ ^ 1 )+( 2/C2_ ^ 2 ). 
the edges of the graph are colored red and black, with 2k\ — £\ black edge 
multiplicity and 2/^2 — £2 red edge multiplicity. Let M(k,£) be the bases of the 
(fc, -£)-sparsity matroid; then (fei, £\, &2, -^-nested sparsity is the intersection of 
the following two matroids, defined by their bases: 

1. Mi = {E' C E\E' e M(k u £ 1 )}, bases in the M(fci, ^-sparsity matroid 
when edge color is disregarded. 

2. M 2 = M(k 2 j2) UKf r!l , bases in the red (fc 2 , ^-sparsity matroid 
padded with full edge multiplicity of the black edges. 

□ 

As a consequence, when < £{ < 2ki, the matroid intersection algorithm 
of Edmonds [4 can be used to solve the Decision (is a graph nested sparse?), 
Extraction ( given an input graph, extract a maximum- sized nested sparse sub- 
graph) and Components (given an input graph, extract its maximal vertex sets 
that span nested tight subgraphs) problems. 

Edmonds' algorithm outputs a maximum-sized set of edges that are inde- 
pendent in both matroids and requires an oracle to test for independence in 
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each matroid. For the oracles, we use the pebble games algorithms of [14], a 
family of algorithms parametrized by two constants k and l\ the (fc,£)-pebble 
game characterizes (&,£) -sparsity. In particular, the (k, £) -pebble game takes 
a graph as input and can be run in two modes: the Decision mode returns 
"yes" if the input graph is (fc, £)-sparse, and the Components mode returns 
the maximal vertex sets that span (/c,£)-tight subgraphs. Algorithm [I] gives a 
more detailed description of how Edmonds' matroid intersection algorithm is 
used to solve problems for nested sparsity. 



Algorithm 1. &2, £2) -nested sparsity: Decision, Extrac- 

tion and Components 

Input: A graph G = (V,E = R{JB) and constants fci, £\, k^ £2, where 

0<£i< 2ki fori = 1,2. 

Method: 

1. Run Edmonds' matroid intersection algorithm CV on G for the 
two matroids M 1 and M 2 , as defined in the proof of Theorem \6.1\ 
When the algorithm performs independence queries on a set of 
edges I C E, 

(a) For the matroid M\, play the (k\,£\) -pebble game on the in- 
put (V,I) in Decision mode, which returns "yes" if it is 
(&i, £1) -sparse 

(b) For the matroid M2, play the (&2, £2) -pebble game on the in- 
put (V,lnR) in Decision mode, which returns "yes" if it is 
(fc 2 , £2) -sparse 

2. Edmonds ' algorithm returns a set ICE that is of maximum size, 
where (V,/) is (fci, £1, &2, £2) -nested sparse. 

3. Output: 

• For Decision, "yes" if 1 — E and no otherwise. 

• For Extraction, (V,/). 

• For Components, play the (fci, £1) -pebble game in Compo- 
nents mode on ( V, I) and output the components returned by 
the pebble game. 



Figure 17: Algorithm for nested sparsity. 

Complexity analysis. Edmonds' algorithm queries the oracles 0(mr 2 ) times, 
where m is the number of elements in the ground set and r is the smaller 
rank of the two matroids. For nested sparsity, on a graph G with n vertices 
and m = 0{n 2 ) edges, both matroids have rank 0{n). Therefore, Edmonds' 
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algorithm requires 0(n 4 ) oracle queries. The pebble game algorithms require 
0(n 2 ) time, resulting in 0(n 6 ) total complexity for Algorithm[l] We note that, 
using the recent matroid intersection algorithm of Harvey [8], a more efficient 
running time of (^(mr^ -1 ) = 0(n u;+1 ), where uo is the matrix multiplication 
exponent, can be obtained for nested sparsity. 

m " we can 



Since (6, 6, 3, 3)-nested sparsity meets the conditions for Theorem 6.1 



apply Algorithm [T] to address the necessary condition for body-and-cad rigidity 
in polynomial time. 



7 Conclusions 

Constraint-based CAD software contains a rich set of geometric constraints. 
Motivated by their applications, we have initiated the study of body-and-cad 
rigidity by identifying a class of constraints amenable to rigidity-theoretical in- 
vestigation and developing their infinitesimal theory. The shape of the rigidity 
matrix naturally led to the study of (6, 6, 3, 3)-nested sparsity, a necessary, but 
not sufficient, condition for body-and-cad rigidity. The polynomial time algo- 
rithm we presented for testing this condition may have practical applications as 
a filter for finding rigid components in a CAD environment, providing informa- 
tive feedback to the user. 

Applications. The results presented can be applied to a larger set of CAD con- 
straints via simple reductions. In particular, it is easy to establish the following 
reductions: 

• Sphere-sphere tangency: Reduces to point-point distance using the 
sphere centers and the sum of the radii. 

• Sphere-plane tangency: Reduces to point-plane distance using the 
sphere center, the plane and the sphere radius. 

• Sphere-line tangency: Reduces to point-line distance using the sphere 
center, the line and the sphere radius. 

• Sphere-point coincidence: Reduces to point-point distance using 
the sphere center, the point and the sphere radius. 

Analogous reductions can be applied when considering cylinders instead of 
spheres by substituting the cylinder's center axis for the sphere's center point. 

Future directions. It remains an open problem to find a combinatorial charac- 
terization for generic body-and-cad rigidity. We anticipate the study of some of 
the constraints introduced here may prove more tractable than classical 3D bar- 
and-joint rigidity. A full combinatorial characterization for angular constraints 
appears in [16J[T3]. However, we observe that finding a complete characteriza- 
tion may require overcoming well-known obstacles such as detecting dependen- 
cies in 3D bar-and-joint, 2D points-and-angles, 2D circles-and-angles, and 2D 
point-line incidence constraint systems. 
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Analogous body-and-cad structures for 2D consist of rigid bodies with pair- 
wise coincidence (point-point, point-line and line-line), angular (line-line) and 
distance (point-point, point-line and line-line) constraints identified between 
points and lines rigidly attached to bodies. The development of the rigidity 
matrix is a straightforward extension of this work. The interesting question, 
which remains future work, is a full combinatorial characterization. 
Acknowledgements. The authors would like to thank Louis Theran for the 



observation and proof of Theorem 6.1 The figures were made using the Solid- 
Works CAD software [3]. 
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